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$8 3e. 近 地 层 空气 动力 学 参数 是 风沙 研究 的 重要 内 容 。 雅 鲁 藏 布 江 (以 下 简称 雅 江 ) 中 游 是 西藏 地 区 风沙 灾害 


最 为 严重 的 区 域 之 一 , 河 漫 浴 和 山 荐 流沙 区 是 雅 江 中 游 地 区 风沙 灾害 的 主要 物 源 区 ,对 其 近 地 表 空气 动力 学 参数 


的 研究 ,有 助 于 对 该 地 区 的 风沙 灾害 进行 评估 和 阔 尘 暴 预 报 预警 。 选 取 河 漫 滩 植被 区 (植被 盖 度 >20%) 人 为 干扰 
区 (植被 盖 度 <20% ) 和 流沙 区 (无 任何 植被 )3 种 下 垫 面 , 利 用 风速 廊 线 法 计算 摩 阻 风速 4) 和 空气 动力 学 粗糙 度 


(zo)。 结 果 表明 :(1) 近 地 层 风速 降幅 植被 区 > 人 为 干扰 区 > 流沙 区 ,植被 区 风速 降幅 约 是 流沙 区 的 2 倍 。(2) 不 同 下 


执 面 近 地 层 风速 廊 线 均 可 以 表达 为 对 数 线性 函数 (w=bln(z)+a) ,但 拟 合 参数 不 同 。(3) 植被 区 、 人 为 干扰 区 和 流沙 区 


平均 值 分 别 为 0.73 ms 、0.59 m-s 和 0.36 ms 2,235 (84) 9173 7.1X 107 m ,3.7X107 m 1814.7x10* mo TES DC uF 
aa 约 为 流沙 区 的 2 倍 和 150 倍 。(4) 植被 区 .人 为 干扰 区 和 流沙 区 的 w 均 随 平 均 风 速 增 大 而 增 大 , 且 增幅 植被 区 > 人 为 
干扰 区 > 流沙 区 ,(5) 植被 区 和 流沙 区 za 对 平均 风速 的 响应 方式 不 同 ,植被 区 的 za 随 风 速 增 大 而 缓慢 增 大 ,流沙 区 的 
zo 在 低 风速 有 减 小 趋势 , 当 风 速 增 大 到 闷 值 (6.5 m:s”) 后 ,za 开始 增加 。 

关键 词 : 雅鲁藏布江 ; 风速 廊 线 ， 摩 阻 风速 ， 空 气动 力学 ;粗糙 度 长 度 


风沙 活动 是 影响 西北 干旱 . 半 干 旱 区 生态 环境 
与 社会 可 持续 发 展 的 主要 因子 之 一 "。 青 藏 高 原 地 
区 同样 存在 风沙 活动 , 因 高 海拔 、 低 气压 的 大 气 环 境 
导致 其 风沙 活动 过 程 与 西北 低 海拔 地 区 明显 不 同 ”。 

空气 动力 学 参数 摩 阻 风速 (w) 和 空气 动力 学 粗 
糙 度 (za) 是 风沙 活动 研究 的 重要 内 容 ”” 。 正 确 计算 
必 和 za, 有 利于 预测 和 量化 输 沙 量 和 粉尘 通 量 ””、 土 


丘 等 11 种 下 垫 面 进行 空气 动力 学 参数 研究 , 发现 沙 
丘 地 表 必 和 aa 最 小 , 且 必 和 az 与 下 垫 面 的 生物 量 以 及 
起 伏 程 度 密切 相关 。 何 玉 斐 等 "认为 ,戈壁 地 表 粗 
PEREN 107-10? m。 周 杰 等 研究 了 策 勒 县 高 风速 
事件 中 沙漠 -绿洲 过 渡 带 4 种 下 垫 面 的 wz, 结果 表 
明 必 和 aa 流动 沙丘 > 半 流 动 沙丘 > 固定 沙丘 > 绿洲 。 
上 述 研 究 表明 ,不 同 下 垫 面 空气 动力 学 性 质 差 异 显 


Epi 7] 7 ,是 评价 防 沙 治 沙 工程 的 重要 指标 ”。 
ww 是 气流 对 地 表 剪 切 力 的 反映 , 当 摩 阻 风 速 超过 临 
界 摩 阻 风速 时 ,地 表 微粒 脱离 地 面 进入 空中 ,发 生 
风蚀 ”“。az 是 地 表 风 速 减 小 到 零 时 的 某 一 高 度 ,是 大 
气 边界 层 水 流 属性 通 量 参数 化 方案 中 常用 的 参数 
之 一 "。az 对 地 表 粗 糙 元 变化 的 啊 应 敏感 ,反映 地 表 
粗糙 元 对 风速 的 减弱 作用 ”。vw aa 受 下 垫 面 性 质 影 
响 显 著 。 茅 宇 壹 等 ”对 不 同 干扰 程度 的 草原 以 及 沙 
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著 。 因 此 ,对 风沙 区 不 同 下 垫 面 近 地 层 气流 动力 学 
参数 的 观测 与 计算 有 利于 认识 不 同 粗 糙 元 与 风 动 
力 之 间 的 能 量 转换 关系 ,进而 评估 区 域 风沙 活动 强 
BE .估算 沙 尘 释放 量 等 。 然 而 ,上 述 研 究 主要 是 针 
对 低 海 拔 的 风沙 区 ,高 海拔 地 区 空气 密度 (0.84 kg， 
m”) 小 于 低 海拔 地 区 (1.23 kg* m?) ;高 海拔 地 区 风沙 
沉积 物 富 含 粉 沙 和 黏土 (5.8%+7.2%)52 , 远 高 于 低 
海拔 地 区 (<1% ) 。 根 据 Shao 等 所 的 起 动 风 速 计算 
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公式 ,高 海拔 地 区 沙 粒 的 起 动 过 程 不 同 于 低 海拔 地 
区 。 但 目前 关于 高 海拔 地 区 的 野外 实测 研究 很 少 ， 
限制 了 我 们 对 高 海拔 地 区 风沙 灾害 形成 机 理 的 
认 知 。 

雅 江 中 游 风沙 活动 强烈 ”, 是 除 中 国 北方 蕊 并 
区 外 的 重要 粉尘 源 区 之 一 汪 。 其 风 旱 同期 的 气候 
和 广 布 的 松散 沉积 物 ”以 及 人 为 破坏 地 表 植 被 与 
土 层 结构 等 ” ,造成 该 区 域 风沙 灾害 严重 ,但 对 该 
区 域 风 沙 问 题 的 研究 相对 薄弱 “。 目 前 , 雅 江 风沙 
研究 主要 集中 在 以 下 几 方 面 :(1) 沉积 物 物 源 。 研 
究 认 为 河谷 底部 ( 河 心 洲 、 河 漫 滩 ) 是 发 生 风 沙 活动 
的 物 源 (2) 风沙 地 貌 。 风 沙 地 貌 类 型 复杂 
多 样 “…" ,但 最 为 典型 的 是 怜 坡 沙丘 。(3) 沙漠 化 。 
该 区 域 土 地 沙漠 化 是 气候 变化 与 人 为 活动 共同 作 
用 的 结果 0.04) 风沙 运动 特征 。 不 同 地 表 风 沙 运 
动 特征 不 同 ””。 但 上 述 研 究 不 能 很 好 阐明 不 同 地 
表 的 沙 全 运动 机 理 , 主 要 原因 是 对 不 同 下 垫 面 的 空 
气动 力学 性 质 的 研究 较 少 。 因 此 ,本 文通 过 野外 实 
测 的 风 廊 线 数据 计算 了 植被 区 、 人 为 干扰 区 和 流沙 
区 3 种 下 垫 面 的 空气 动力 学 参数 ,为 区 域 输 沙 量 计 
算 ,粉尘 通 量 及 区 域 风沙 灾害 评估 提供 参考 。 


1 研究 区 概况 


研究 区 位 于 雅鲁藏布江 ( 雅 江 ) 中 游 山 南 段 ， 


西 起 贡 嘎 县 , 东 至 乃 东 区 (90°30'~92°40'E、29°10'~ 
29°30'N) ,海拔 高 于 3550 m, 上 且 自 西向 东 逐 渐 降 低 
(图 1) , 北 侧 为 念 青 唐古拉 山 , 南 侧 为 喜马拉雅 山 。 
研究 区 年 平均 气温 为 8.7 % ,年 平均 降水 量 为 378 
mm ,年 内 变化 曲线 呈 单 峰 型 , 旱 ( 冬 .春季 ) LR CELL) 
季 分 明 ;年 平均 相对 湿度 为 429 ,夏季 较 湿 .冬季 和 较 
干 ;年 平均 风速 为 2.6 ms, KEER, RIKEN; AK 
向 复杂 , 受 地 形 影 响 显著 , 主 风向 为 西风 。 

雅 江 中 游 水 系 兴 状 发 育 , 受 径流 季节 性 影响 ， 
冬 春季 水 位 低 , 河 心 洲 及 河 漫 滩 大 面积 裸露 7， 
表层 干燥 松散 的 沉积 物 在 强 气 流 作用 下 向 北面 山 
苞 甚 至 山坡 输 运 ,风沙 活动 频繁 '"”。 有 关 人 研究 表 
HH ,该 地 区 的 风沙 灾害 可 能 始 于 以 青 牺 农 耕 为 代表 
的 历史 时 期 ,并 随 着 人 类 开发 强度 的 加 大 ,风沙 灾害 
日 趋 严 重 ;近年 来 ,风沙 活动 频次 有 增多 的 趋势 ”。 


2 研究 方法 


2.1 野外 测量 

人 研究 区 冬 .春季 风沙 活动 频繁 , 尤 以 正午 之 后 
更 为 严重 , 故 选 择 2021 年 1 月 26 日 14:35 一 18:04 和 
1 月 27 日 12:35 一 1$:34 分 别 位 于 山南 市 贡 咕 县 和 森 布 
日 (图 1]b) 和 昌 果 附近 的 河 漫 滩 处 (图 1e ) ,利用 二 维 
超声 波 风速 仪 (风速 量程 分 别 为 0~30 m-s”, DAR 


29?20'N 


mu 河 心 洲 林地 
mm 流动 沙 地 o 城镇 
e 测 点 


um 河流 
m 河 心 洲 
m 半 固 定 沙 地 固定 沙 地 


29?10'N 


(b) 森 布 日 


为 0.01 ms ;风向 量程 为 0~359° ;分 辩 率 为 1°; 数 
91?20'E 91?40'E 
Z 
S 
& 
Z 
S 
PA E 


图 1 研究 区 地 理 位 置 及 各 站 点 地 表 状 况 


Fig. 1 Location of study area and sites and land surface properties in Semburi and Changguo, respectively 
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据 记 录 频 率 为 1 min) 对 不 同 高 度 风速 进行 观测 。 
观测 期 间 ,两 地 均 形 成 明显 的 风沙 流 。 

条 布 日 地 表 有 植被 覆盖 ,以 芦苇 (Phiragrmiies 
australias ) JJ E ,植被 高 度 为 20~40 cm ,植被 覆盖 度 
为 10%~30%( 图 lb)。 昌 果 地 表 为 松散 的 风沙 沉积 
物 , 无 植被 覆盖 (图 1c)。2 个 测 点 地 表 起 伏 程度 均 
较 小 。 在 森 布 日 选择 植被 区 (S1 ,植被 盖 度 >20% ) 和 
人 为 干扰 区 (S2 ,植被 盖 度 <20%)2 种 下 垫 面 进行 测 
量 ; 设 有 5 层 高 度 ,风速 传 感 融 安 装 高 度 分 别 为 8 
cm,20 cm,40 cm,90 cm 和 180 cm。 虽 然 , 昌 果 下 垫 
面 性 质 相 近 ,但 为 了 对 比 ,同样 测量 了 2 组 风 廓 线 
《C1,C2; 均 无 植被 ) , 且 风 速 传 感 顺 安装 高 度 分 别 为 
5 cm, 10 cm 、 30 cm, 80 cm,170 cm。 两 地 参考 风速 
的 风速 传感器 架 均 设 在 距 地 面 200 em ir FEM, fk 
布 日 和 昌 果 分 别 记录 210 条 和 180 条 数据 。 

2.2 数据 分 析 

高 风速 条 件 下 大 气 层 结 呈 中 性 ” ,中 性 层 结 条 

件 下 近 地 层 风 廊 线 方程 可 以 用 下 式 ” 表 示 : 

ufu. = In(z/zo)/k (1) 
式 中 :表示 高 度 z 处 的 风速 (ms ') su ZI HL EE 
阻 风速 (ms ');z 表 示 高 度 (m);% 表 示 空 气动 力学 
粗糙 长 度 (m);k 是 冯 : 卡 曼 常 数 (k=0.4)。 

风速 与 高 度 的 关系 可 以 按照 下 式 " 利用 最 小 二 
乘法 进行 回归 拟 合 。 同 时 ,本 文 引入 判定 系数 (RR) 
和 均 方 根 误差 (RMSE ) 以 评价 拟 合 优 度 。 

u,=bln(z)+a (2) 


森 布 日 (1 月 26 日 14:35 一 18:04) 


(a) . 
20 boe a x uad T ^N AE 
7 n fa! Put Te ie Sa 


风速 /ms 


时 间 /min 


图 2 观测 其 


式 中 ;a,b 为 拟 合 系数 。 
ww 和 zw 可 分 别 用 下 式 计算 : 
u. = kb (3) 
zo = exp(— a/b) (4) 
在 用 风 廓 线 方程 计算 参数 时 ,为 了 控制 数据 质 
量 , 选 择 2m 高 度 风速 >4m.s- H.JXUS] 2g 2009-290? 
范围 内 的 数据 (图 2 中 阴影 区 域 )。 其 中 , 森 布 日 有 
158 > 4 m- s” 的 风速 记录 (平均 风速 为 7.20+2.17 
m's, EKUB W 24821?) , BRA 1734824 m s" 
的 风速 记录 (平均 风速 为 7.98+1.56 ms ', 主 风向 为 
240°+14°)。 
为 了 对 比 不 同 下 垫 面 的 风速 廊 线 ,将 不 同 高 度 
的 风速 和 最 高 处 的 风速 按照 下 式 进行 归 一 化 处 理 : 
u; —-uju, (5) 
式 中 :i 表示 从 低 到 高 第 i 层 (i=1,2,3,4,5);w 是 第 i 
E USE (m «8 ;w' 是 第 i 层 风 速 归 一 化 之 后 的 无 量 
纲 风 速 ;wW 是 参考 风速 (ms )。 


3 结果 与 分 析 


3.4 风速 廓 线 特 征 

植被 区 和 流沙 区 风速 随 高 度 降低 而 降低 ,但 风 
速 降幅 差异 显著 ( 表 1)。 各 站 点 风速 降幅 为 :植被 
区 > 人 为 干扰 区 > 流沙 区 (C1、C2) , 413] 2 : 5296 
21% 40%+19% 18%+12% .18%+13%, 

2 种 下 垫 面 的 风速 廓 线 均 符 合 对 数 线性 函数 关 


昌 果 (1 月 27 日 12:35 一 15:34) 


1 
30 60 90 120 150 180 


风速 /ms 


时 间 /min 


注 :阴影 区 域 为 论文 分 析 数 据 。 
间 风 速 和 风向 


Fig.2 Wind speed and direction during field experiments 
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表 1 和 森 布 日 和 昌 果 不 同 高 度 的 平均 风速 及 风速 降幅 


Tab.1 Averaged five-layer wind speeds and its reduction amplitude at Semburi and Changguo 


高 度 /cm ”SI1 风 速 /(m*s") 降幅 /% ”5S2 风 速 /(m:s"') ”降幅 /% | 高度/cm C1 风速 /(ms"') ”降幅 /% ”C2 风速 /(m:s"') ”降幅 /% 
180 6.43+2.19 31 6.36+2.27 30 170 8.12+1.61 7 7.77+1.67 6 
90 4.44+1.97 36 4.45+2.05 32 80 7.59+1.44 9 7.30+1.62 16 
40 2.86+1.47 45 3.02+1.35 29 30 6.88+1.27 12 6.12+1.51 5 
20 1.56+0.98 58 2.15+1.18 31 10 6.07+1.08 21 5.83+1.27 21 
8 0.65x0.37 89 1.49+0.79 TT 5 4.80+0.78 40 4.58+1.07 41 
注 :S1( 植 被 区 ,植被 盖 度 > 20% ) 为 森 布 日 1 点 ;S2( 人 为 干扰 区 ,植被 盖 度 <20%) 为 森 布 日 2 点 ;C1 为 昌 果 1 点 ;C2 为 昌 果 2 点 。 下 同 。 


系 (w=bln(z) + a, R>0.93, RMSE 为 0.41-0.52, 图 
3a)。 下 垫 面 性 质 不 同 拟 合 函 数 的 参数 明显 不 同 。 
风速 廊 线 的 斜率 越 大 ,风速 随 高 度 降低 的 幅度 越 
大 。 植 被 区 相对 人 为 干扰 区 植被 高 度 高 .覆盖 度 大 ， 
BEST 点 和 斜率 (1.83) 大 于 S2 点 斜率 (1.47)。 昌 果 地 表 
均一 , 均 为 流沙 , 故 昌 果 两 点 斜率 大 小 相近 (0.91、 
0.87)。 归 一 化 的 风速 廓 线 同样 表现 出 相似 的 规律 
CE] 3b) :植被 高 度 越 高 .覆盖 度 越 大 , 则 斜率 越 大 ， 
S1 点 和 S2 点 和 斜率 分 别 为 0.28 .0.23。C1 点 和 C2 点 
斜率 均 为 0.11 ,说 明 本 研究 的 结果 具有 可 重复 性 , 风 
速 廊 线 系数 的 差异 是 由 下 热 面 属性 不 同 造成 的 。 
3.2 空气 动力 学 参数 特征 

下 热 面 性 质 不 同 , 必 差 异 显著 。 以 在 植被 区 > 人 
为 干扰 区 > 流沙 区 (图 4a, 图 4b)。 植 被 区 已 约 为 人 
为 干扰 区 的 1.2 倍 , 约 是 流沙 区 的 2.1 倍 。 植 被 区 六 
的 分 布 范 围 广 , 呈 单 峰 分 布 ,存在 粗 尾 (图 4a)。S1 
点 的 必 最 小 值 .平均 值 和 最 大 值 分 别 是 0.20 m-s, 
0.73 m*s $1.76 m-s (图 4b)。S2 点 的 必 最 小 值 、 
平均 值 和 最 大 值 分 别 是 0.19 ms 0.59 m s^ ff 
1.54 m-s% FLW DX u AYANA V FE] FEER AIPE , u- 


- (a) 实测 风速 廓 线 


Inz/m 


0.5 2.5 4.5 6.5 8.5 


u/(m:s?) 


分 布 集中 。C1 ELT] uds] ME .平均 值 和 最 大 值 分 别 
是 0.19 m-s,0.35 m-s” 和 0.69 m-s C2 IU ufi 
小 值 .平均 值 和 最 大 值 分 别 是 0.19 m.s 0.36 m-s” 
和 0.67 ms 。 

下 垫 面 性 质 不 同 ,am 差 异 亦 显著 (图 4c, 图 4d)。 
2 在 植被 区 > 人 为 干扰 区 > 流沙 区 , 且 植 被 区 za 约 是 
人 为 干扰 区 的 2 倍 , 约 是 流沙 区 的 150 倍 。 植 被 区 zo 
的 变化 区 间 大 ,概率 密度 曲线 存在 细 尾 (图 4c)。S1 
HI z "E EX [8] JJ :0.5x107-13.1x10? m( 图 4d) ,平均 
值 为 7.1x10? m; 人 为 干扰 区 (S2)zo 变 化 区 间 为 :0.1x 
107-10.4x10? m, 平 均值 为 3.7x10” mo Tat Ub EX zo HY 
变化 区 间 比 植被 区 小 ,zo 分 布 较 植 被 区 集中 ,近似 单 
峰 分 布 。C1 的 z 变 化 区 间 为 :1.1x10“~39.3x10” m, 
平均 值 为 3.2x107 m; C2 的 % 变 化 区 间 为 :1.2x10“~ 
47.8x10* m, 平 均值 为 3.5x10™ m. 


4 讨论 


41 空气 动力 学 参数 与 下 垫 面 性 质 的 关系 
空气 动力 学 参数 受 地 表 属性 和 气流 状态 共同 
影响 ,其 变化 是 粗糙 元 与 气流 相互 作用 的 结果 ”。 


cp 


(b) 归 一 化 风速 廓 线 


图 3 KERR 
Fig.3 Wind profile 
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注 :图 a, 图 b 为 森 布 日 和 昌 果 uw; 图 6, 图 d 为 森 布 日 和 昌 果 z。 
图 4 摩 阻 风 速 (w) 和 空气 动力 学 粗糙 度 (z) 概 率 密 度 


Fig.4 Probability density curves of u» and Zo 


为 了 阐明 下 垫 面 对 aa 和 六 的 影响 ,本文 计 算 参 数 时 
不 仅 严格 限制 了 风速 和 风向 (图 2) ,还 对 风速 进行 
了 归 一 化 处 理 (图 3b)。 结 果 表明 ,植被 区 和 流沙 区 
的 空气 动力 学 参数 差异 显著 (P<0.05), 反 映 了 下 热 
面 性 质 对 气流 状态 影响 显著 。 植 被 区 和 流沙 区 地 
表 主 要 差异 在 于 植被 盖 度 和 植被 高 度 ,植被 区 植被 
TREAN 1096-3096 ,植被 高 度 约 为 20~40 cm ,流沙 
区 植被 盖 度 为 0( 图 16) ;其 实质 反映 了 地 表 粗 糙 元 
的 不 同 。Lettau 等 2 认为 影响 地 表 空 气动 力学 参数 
大 小 的 原因 很 复杂 ,但 最 关键 的 决定 因素 是 下 垫 面 
的 粗糙 元 。 瑟 等 呈 得 出 相似 结论 ,认为 地 表 粗 糙 元 
是 影响 空气 动力 学 参数 的 关键 因子 ,地 表 粗 糙 元 越 
大 ,uw 和 如 愈 大 。 植 被 的 存在 增加 了 地 表 粗 烽 元 , 增 
强 了 近 地 层 气流 的 雷诺 剪 切 引力 且 造 成 零 平 面 高 
度 上 升 到 为, 故 植被 区 的 性 和 za 大 于 流沙 区 。 如 S1 
植被 盖 度 和 高 度 明 显 大 于 S$2 ,导致 S1 HY zoze S2 的 
1.9 倍 。C1 和 C2 为 无 植被 ,所 以 a 平均 值 变 化 不 明 
显 ,分 别 为 3.2x10“*m 和 3.5x10“*m。 植 被 区 的 顶层 
和 底层 风速 降幅 可 达 89% ,而 流沙 区 仅 为 40%。 这 
种 差异 导致 植被 区 的 明显 大 于 流沙 区 (SI1 和 C1 的 
忆 平 均值 分 别 为 0.73 mes” 410.35 mes) SEIREN 
认为 造成 这 种 差异 的 主要 原因 是 柔性 植被 可 以 通 
过 改变 自身 的 结构 和 形态 更 多 地 消耗 近 地 表 气流 
的 能 量 ,减少 了 直接 到 达 地 表 的 能 量 ,从 而 增 大 了 
近 地 表 风速 的 降幅 ,使 得 近 地 表 风 速 快 速 降 至 0, 导 
致 必 和 z 增 大 。 


4.2 空气 动力 学 参数 与 平均 风速 的 关系 

4.2.1 摩 阻 风速 与 平均 风速 的 关系 ， 植 被 区 和 流沙 
区 的 与 2 m 高 度 平均 风速 呈 显著 的 线性 正 相 关 关 
系 (P<0.05) ,这 与 前 人 的 研究 结果 相似 ” ,但 2 种 
下 肆 面 的 心 随 风 速 增 大 程度 不 同 (图 5)。 植 被 区 六 
随 风速 增 大 程度 是 流沙 区 的 2 倍 (S2 和 C2 斜率 分 别 
为 0.12 和 0.06)。 此 外 ,S1 增 加 幅度 略 大 于 S2( 和 斜率 
分 别 为 0.13 和 0.12) , C1 与 C2 地 表 均 是 由 平坦 的 流 
沙 组 成 ,下 垫 面 相似 ,两 地 uw: 随 风 速 增 大 的 程度 相 
近 。 据 此 判断 , 必 与 风速 的 增 大 程度 与 地 表 类 型 密 
切 相关 ,植被 覆盖 的 地 表 w 随 风速 增 大 程度 比 流沙 


ums) 


4 6 8 10 12 


u/(ms?) 


图 5 摩 阻 风速 (x,) 与 平均 风速 (z) 的 关系 
Fig.5 Relationship between threshold friction velocity (u-) 


and averaged wind speed (u) 
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地 表 大 ; 近 地 层 气流 滑 流 强度 增强 是 地 表 摩 阻 风速 
增 大 的 主要 原因 。 在 相同 地 表 , 近 地 层 消 流 强度 随 
着 风速 增 大 而 增 大 ,导致 地 表 摩 阻 风速 增 大 ,表现 
出 摩 阻 风 速 随 风速 增 大 而 增 大 。 除 此 之 外 ,笔者 还 
观察 到 植被 区 摩 阻 风速 随 风速 的 增 大 程度 比 流沙 
区 大 ;其 原因 是 粗糙 度 大 的 植被 区 增强 了 近 地 层 气 
流 的 潮流 强度 , 且 随 风速 增 大 而 增强 。 
4.2.2 空气 动力 学 粗糙 度 与 平均 风速 的 关系 ， 植 被 
区 和 流沙 区 za 与 2 m 高 度 平均 风速 的 关系 明显 不 同 
(图 6)。 植 被 区 的 随 风速 增 大 而 缓慢 增 大 (图 
6a) 。 流 沙 区 的 aa 随 风速 增 大 先 缓慢 减 小 , 当 风 速 
JE KI] (EL (6.5 m «s^ ') 后 ,zw 迅速 增 大 (图 6b)。 张 
强 等 5 通过 研究 不 同 植被 下 垫 面 的 空 力 动 力学 参 
数 时 发 现 , 近 地 层 风速 会 明显 改变 柔性 粗糙 元 的 结 
构 和 形态 ,从 而 导致 % 在 低 风 速 时 较 大 ,高 风速 时 较 
小 。 杨 兴 华 等 ?在 塔克拉玛干 沙 漠 过 渡 带 研究 了 
无 风沙 运动 的 床 面 表 明 ,大 和 气 层 结 呈 中 性 时 ,za 与 2 
m 高 度 处 风速 旺 显 著 的 负 指数 关系 。Zhang 等 ”在 
腾 格 里 沙漠 东南 缘 观 测 风沙 活动 条 件 下 , 草 方 格 履 
盖 的 地 表 和 沙 质地 表 的 空气 动力 学 参数 时 发 现 , 低 
风速 (<6 m:'s-) 时 草 方 格 覆 盖 的 地 表 a 随 风速 迅速 
增 大 ,高 风速 (>6 m.s) 时 缓慢 增 大 ;但 沙 质 地 表 a 


(6.5 m*s ') 指 示 了 起 动 风速 ,这 与 计算 的 临界 起 沙 
风 风 速 相似 (流沙 地 表 平 均 粒 径 为 292.81 pm。 沙 
粒 密度 为 2650 kg*m”。 观 测 期 间 平均 气压 为 655.7 
hPa, 平 均 气 温 为 10.1 C。 根 据 气 体 状 态 方程 计算 
出 空气 密度 约 为 0.81 kg.m3?。 据 此 ,利用 起 沙 风 公 
式 引 计算 流沙 地 表 2 m 高度 处 临界 起 沙 风 风速 为 
6.6m.s-)。 
4.3 海拔 对 风速 廊 线 的 影响 

为 了 阐明 海拔 与 近 地 层 风速 的 关系 ,对 不 同 区 
域 流沙 地 表 的 野外 实测 资料 按照 公式 (5) 做 归 一 化 
处 理 ( 图 7)。 研 究 区 流沙 地 表 ( 昌 果 ) 与 塔 中 2 X 
煌 字库 布 齐 风电 场 中 ` 巴 丹 吉 林 拐 子 湖 六 .中卫 风 
沙 观 测 场 3 的 归 一 化 风速 廊 线 均 符合 对 数 线形 规 
律 (w'=cln(z)+d;c,d 为 标准 化 风速 廓 线 拟 合 系数 ， 
R°>0.98 , RMSE 为 0.03~0.06) (KI 7a, 3€ 2) ;参数 c 随 
海拔 增加 呈 指 数 递减 ,但 & 呈 指数 递增 ,揭示 了 近 地 
层 风速 随 海拔 高 度 变化 的 差异 (图 7b .图 7c);c 值 越 
大 (d 值 越 小 ) ,风速 减 小 越 快 ,反之 则 亦 然 。 这 表明 
相对 于 低 海 拔 地 区 ,高 海拔 地 区 空气 密度 小 ,黏度 
小 ,气流 内 摩擦 力 小 ,气流 能 量 向 下 传递 时 被 消耗 
的 少 , 根 据 雷 诺 数 公式 判断 高 海拔 地 区 更 易于 发 生 
满 流 。 说 明 高 海拔 地 区 气流 随 高 度 增加 ( 减 小 ) 程 度 


在 低 风速 (<8 m's”) 时 随 风速 增 大 而 缓慢 增 大 ,在 高 
风速 (>8 ms ) 时 aa 迅速 增 大 。 高 咏 晴 等 ”通过 风 
洞 实验 研究 了 风沙 流 和 净 风 场 中 空气 动力 学 参数 
发 现 , 风 沙 流 中 z 大 于 净 风 场 , 且 甩 在 净 风 场 中 随 着 
风速 增 大 而 减 小 ,但 在 风沙 流 中 z 随 着 风速 的 增 大 
而 增 大 。 王 染 等 ”的 研究 也 证 明了 这 一 点 。 综 上 
所 述 , 笔 者 认为 :(1) 植被 区 za 是 风沙 流 和 植被 综合 
作用 的 结果 。(2) 流沙 区 有 无 风沙 流 对 计算 的 粗粮 
度 值 影响 显著 。 无 风沙 流 时 计算 的 za 可 能 会 俩 大 ， 
形成 风沙 流 时 ,zw 随 风速 呈 指 数 增加 。 而 临界 值 


(a) Rt H @ SI A S2 


e 
E 6 | zj-0.06e29» RI-0.54 RMSE-0.009 
E z-0.04e^0« R?=0.42 RMSE=0.005 
^ ^ 


u/(m's?) 


要 比 低 海拔 地 区 慢 , 能 量 消耗 较 少 ,从 而 导致 高 海拔 
地 区 风沙 沉积 物 比 低 海拔 地 区 运动 的 更 高 .更 远 ”。 

综 上 所 述 , 森 布 日 和 上 昌 果 的 性 和 aa 差异 显著 ,与 
平均 风速 的 关系 明显 不 同 。 尽 管 在 高 风速 条 件 下 ， 
森 布 日 和 昌 果 地 表 都 有 可 能 发 生 风 蚀 , 但 通过 对 和 森 
布 日 和 昌 果 凡 和 aa 的 分 析 , 昌 果 地 表 发 生 风 蚀 的 湾 
力 更 大 。 结 合 研究 区 的 地 表 类 型 时 空 变 化 特征 ， 
A .春季 雅 江 中 游 类 似 昌 果 地 表 环 境 的 区 域 , 河 漫 
滩 河 心 洲 大 面积 裸露 ” ,建议 风沙 灾害 防治 重点 
放 在 河 心 洲 和 河 漫 滩 。 
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Fig. 6 Relationship between aerodynamic roughness length (zo) and averaged wind speed (u) 
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图 7 归 一 化 风速 廓 线 ,风速 廓 线 参数 (c, d) Eid Ch WIS S 


Fig.7 Normalized wind speed profiles of naturally flat bare sand surfaces at high and low 


altitudes areas and coefficient c and d versus altitude 


32 归 一 化 风速 廊 线 参数 (zw=cln(z)+d) 
Tab.2 Normalized wind speed profile 


parameters for each station 


地 点 海拔 /m c d R^ — RMSE 
pe" 1099 023 0.88 0.98 0.04 
ARWR 1139 0.16 0.92 0.99 0.08 
ti FAREA” 1200 0.13 0.93 0.99 0.03 
巴 丹 吉林 拐子 湖 W 1390 0.12 0.94 0.98 0.06 
中 卫 风 沙 观测 场 和 W1650 0.1 0.94 0.98 0.05 
BR 3650 0.1 0.96 0.99 0.8 

5 结论 


(1) 植被 区 近 地 层 风速 的 降低 幅度 大 于 流沙 
区 ,体现 了 植被 对 近 地 层 风速 的 减弱 作用 ,从 而 减 
弱 土 壤 风 蚀 。 

(2) HERREY us zo KF WUD D ,植被 区 的 uz 
分 别 约 是 流沙 区 的 2 倍 `15$0 倍 。 植 被 区 和 流沙 区 的 
凡 平 均值 分 别 为 0.73 m-s” 和 0.36m.s ,对 应 am 平均 
值 分 别 在 10” 10" 量 级 。 据 此 ,风沙 灾害 防治 重点 
应 该 在 河 心 洲 和 河 漫 滩 。 

(3) 植被 区 与 流沙 区 的 摩 阻 风速 均 与 平均 风速 
呈 显 著 的 线形 关系 , 随 平均 风速 增 大 而 增 大 ,但 植 
被 区 增 大 的 程度 大 于 流沙 区 。 

(4) 植被 区 与 流沙 区 的 空气 动力 学 粗糙 度 对 平 
均 风 速 的 响应 方式 不 同 , 与 下 垫 面 粗糙 元 性 质 以 及 
是 否 形成 风沙 流 密切 相关 。 


参考 文献 (References ) : 


[] ERE, 严 平 , RE. 中 国 北方 风沙 活动 的 驱动 力 分 析 [. 第 四 


[2] 


[3] 


[4| 


[5] 


[6] 


" 


[8] 


[9] 


[10] 


纪 人 研究 , 2001, 21(1): 41-47. [Shi Peijun, Yan Ping, Yuan Yi. The 
driving factors of the blown-sand activity in Northern China[J]. Qua- 
ternary Sciences, 2001, 21(1): 41-47. ] 

Dong Z B, Hu G Y, Qian G Q, et al. High-altitude aeolian research 
on the Tibetan Plateau[J]. Reviews of Geophysics, 2017, 55(4): 
864—901. 

Bagnold R A. The Physies of Blown Sand and Desert Dunes[M]. 
New York: Methuen, 1941: 85-95. 

Dupont S, Rajot J L, Labiadh M, et al. Aerodynamic parameters 
over an eroding bare surface : Reconciliation of the law of the wall 
and eddy covariance determinations[J]. Journal of Geophysical Re- 
search: Atmospheres, 2018, 123(9): 4490-4508. 

Stull R B. An Introduction to Boundary Layer Meteorology[M]. 
New York: Kluwer Academic Publishers, 1988: 378—381. 

杨 兴 华 , 康永 德 , 周 成 龙 , 等 . 塔克拉玛干 沙漠 土壤 粒度 分 布 特 
征 及 其 对 粉尘 释放 的 影响 四 . 农业 工程 学 报 , 2020, 36(5): 167- 
174. [Yang Xinghua, Kang Yongde, Zhou Chenglong, et al. Char- 


acteristics of soil particle size distribution and its effect on dust 
emission in Taklimakan Desert[J]. Transaction of the Chinese Soci- 
ety of Agricultural Engineering, 2020, 36(5): 167-174. | 

Dong Z B, Liu X P, Wang H T. The aerodynamic roughness with a 
blowing sand boundary layer (BSBL): A redefinition of the owen ef- 
fect[J]. Geophysical Research Letters, 2003, 30(2): 1047. 

RIE. 风沙 地 貌 与 治 沙 工程 学 [MI. 北京 : 科学 出 版 社 , 2003: 10- 
20. [Wu Zheng. Geomorphology of Wind Drift Sand and Their Con- 
trolled Engineering[M]. Beijing: Science Press, 2003: 10-20. | 

茅 宇 豪 , 刘 树 华 , 李 婧 . AS In] P Rar zs 6207045 URS IESU]. 
气象 学 报 , 2006, 64(3): 325- 334. [Mao Yuhao, Liu Shuhua, Li 


Jing. Study of aerodynamic parameters for different underling sur- 
faces[J]. Acta Meteorologica Sinica, 2006, 64(3): 325-334. ] 

fu] EAE, 张 宏 升 , 刘 明 星 , 等 . 戈壁 下 垫 面 空气 动力 学 参数 确定 
的 再 研究 加 . 北京 大 学 学 报 (自然 科学 版 ), 2009, 45(3): 439- 
443. [He Yufei, Zhang Hongsheng, Liu Mingxing, et al. A study on 


202210.00036v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


[111] 


[12] 


[13] 


[14] 


[15] 


[16] 


[17] 


[18] 


[19] 


[20] 


determining aerodynamic parameters over Gobi Desert surface[J]. 
Acta Scientiarum Naturalium Universitatis Pekinensis, 2009, 45 
(3): 439-443. ] 

周 杰 , 雷 加 强 , 郑 则 浩 , 等 . 塔里木 盆地 西南 缘 不 同 下 垫 面 的 摩 
阻 风速 和 空气 动力 学 粗糙 度 [ 中 .水 土 保持 通报 , 2016, 36(1): 1- 
5, 10. [Zhou Jie, Lei Jiaqiang, Zheng Zehao, et al. Friction wind 


velocity and aerodynamic roughness over different underlying sur- 
face types in southwestern Tarim Basin[J]. Bulletin of Soil and Wa- 
ter Conservation, 2016, 36(1): 1-5, 10. ] 

IKAk, 马 月 飞 , 兽 林 , 等 . 基于 沉积 物理 化 性 质 的 雅鲁藏布江 中 
JE 2b y ST 2E]. 中 国 沙漠 , 2021, 41(3): 92-100. [Zhang 


Yan, Ma Pengfei, Zeng Lin, et al. Study on silt and clay prove- 


nance in the Yarlung Zangbo River middle reaches using sedi- 
ment physicochemical characteristics[J]. Journal of Desert Re- 
search, 2021, 41(3): 92-100. ] 

Shao Y P, Lu H. A simple expression for wind erosion threshold 
friction[J]. Journal of Geophysical Research Atmospheres, 2000, 
105(D17): 22437-22443. 

李 森 , 杨 萍 , 董 玉 祥 , 等 . 西藏 土地 沙漠 化 及 其 防治 [M]. 北京 : 科 
学 出 版 社 , 2010: 90-98. [Li Sen, Yang Ping, Dong Yuxiang, et al. 


Land Desertification and Its Control in Tibet[M]. Beijing: Science 
Press, 2010: 90-98. | 

Han Y X, Fang X M, Kang S C, et al. Shifts of dust source regions 
over central Asia and the Tibetan Plateau: Connections with the 
Arctic oscillation and the westerly jet[J]. Atmospheric Envireon- 
ment, 2008, 42(10): 2358-2368. 

HE €, IKEI, 论 珠 群 培 , 等 . 雅鲁藏布江 曲 水 - 泽 当 段 风沙 
活动 动力 条 件 分 析 与 风沙 灾害 防治 建议 四 . 中 国 沙漠 , 2021, 
41(1): 10-18. [Ma Pengfei, Zhang Zhengcai, Lunzhu Qunpei, et al. 


Analysis on the sand transport wind power conditions and sugges- 
tions on the sand disaster preventions in the middle reaches of 
Yarlung Zangbo River, China[J|. Journal of Desert Research, 
2021, 41(1): 10-18. ] 

Zl €, 论 珠 群 培 , IKA, 等 . 雅鲁藏布江 中 游 江 心 洲 ` 河 漫 滩 面 
积 及 其 指示 的 沙 源 特征 加 . 中 国 沙漠 , 2021, 41(3): 25-33. [Ma 
Pengfei, Lunzhu Qunpei, Zhang Yan, et al. Sand supplement char- 


acteristics indicated by the area of river inland, flood plain in the 
riverbed of the Yarlung Zangbo River[J]. Journal of Desert Re- 
search, 2021, 41(3): 25-33. ] 

杨 军 怀 , 夏 敦 胜 , 高 福元 , 等 . 雅鲁藏布江 流域 风 成 沉积 研究 进 
IJ]. 地 球 科 学 进展 , 2020, 35(8): 863-877. [Yang Junhuai, Xia 


Dunsheng, Gao Fuyuan, et al. Aeolian deposits in the Yarlung 
Zangbo River Basin, southern Tibetan Plateau: A brief review[J]. 
Advances in Earth Science, 2020, 35(8): 863-877. ] 

TE, 魏 兴 琥 , 董 玉 祥 , 等 . 西藏 沙漠 化 研究 进展 与 未 来 防 沙 治 
vb LR]. 中 国 科 学 院 院 刊 , 2020, 35(6): 699-708. [Yang Ping, 


Wei Xinghu, Dong Yuxiang, et al. Progress on study desertifica- 


tion research and future combating idea in Tibet, China[J]. Bulle- 
tin of Chinese Academy of Sciences, 2020, 35(6): 699—708. ] 
落 桑 曲 加 , IKA, 马 月 飞 , 等 . 雅鲁藏布江 中 游 不 同 地 表 输 沙 量 


i21] 


[22] 


[23] 


[24] 


[25] 


[26] 


[27] 


[28] 


[29] 


[30] 


[31] 


TEAE[J]. 中 国 沙 漠 , 2022, 422): 6-13. [Luosang Qujia, Zhang 


Yan, Ma Pengfei, et al, Study on the sand transport quantities on 


the different landscapes in the middle area of Yarlung Zangbo Riv- 
er[J]. Journal of Desert Research, 2022, 42(2): 6-13. | 

Lettau H. Note on aerodynamic roughness-parameter estimation on 
the basis of roughness element description[J]. Journal of Applied 
Meteorology, 1969, 8(5): 828-832. 

Li X B, Bo T L. Statistics and spectra of turbulence under different 
roughness in the near- neutral atmospheric surface layer[J]. Earth 
Surface Processes and Landforms, 2019, 44(7): 1460-1470. 
Mohammadreza K, David T, Gerard M. A review on turbulent flow 
over rough surfaces: Fundamentals and theories[J]. International 
Journal of Thermofluids, 2021, 10(1): 1-34. 

John F. Turbulence in plant canopies[J]. Annual Review of Fluid 
Mechanics, 2000, 32: 519-571. 

张强 , 曾 剑 , 姚 桐 . 植被 下 垫 面 近 地 层 大 气动 力 状态 与 动力 学 粗 
焰 长 度 的 相互 作用 及 其 参数 化 关系 叮 . 科学 通报 , 2012, 57(8): 
6477-655. [Zhang Qiang, Zeng Jian, Yao Tong. Interaction of aero- 


dynamic roughness length and windflow conditions and its parame- 
terization over vegetation surface[J|. Chinese Science Bulletin, 
2012, 57(8): 647-655. ] 

程 宏 , 分 学 勇 , 张 春 来 . 摩 阻 风速 与 平均 风速 的 转化 关系 研究 
[J]. 水土 保持 研究 , 2007, 61(2): 133-134, 138. [Chen Hong, Zou 


Xueyong, Zhang Chunlai. A study on the relationship between fric- 
tion wind velocity and mean wind velocity[J]. Research of Soil and 
Water Conservation, 2007, 61(2): 133-134, 138. | 
杨 兴 华 , 37] KKH, 张 瑞 军 , 等 . FE DERI E VP EUR 
带 空气 动力 学 粗粮 度 分析 町 . 干旱 区 研究 , 2012, 29(3): 524- 
528. [Yang Xinghua, Ali Mamtimin, Zhang Ruijun, et al. Analysis 


RE 


on aerodynamic roughness length in a desert transitional zone of 
the Taklimakan desert[J]. Arid Zone Research, 2012, 29(3): 524- 
528. ] 

Zhang C L, Lia Q, Zhou N, et al. Field observations of wind pro- 
files and sand fluxes above the windward slope of a sand dune be- 
fore and after the establishment of semi-buried straw checkerboard 
barriers[J]. Aeolian Research, 2016, 20: 59-70. 

高 咏 睛 , 亢 力 强 , 张 军 杰 , 等 . 风沙 流 和 净 风 场 中 瞬时 水 平 风速 
廓 线 特征 比较 册 . 中 国 沙漠 , 2017, 37(1): 48-56. [Gao Yongq- 


ing, Kang Liqinag, Zhang Junjie, et al. Comparison of instanta- 


neous horizontal wind velocity profile in aeolian sand transport 
and clean air[J]. Journal of Desert Research, 2017(37): 48-56. | 

EAR, FEF, 雷 加 强 , 等 . 不 同 下 垫 面 沙 通 量 估算 一 一 以 策 勤 
Ax Ut - Vb WE SEBETE 2 DILE]. 干旱 区 地 理 , 2017, 40(3): 533-540. 


[Wang Cui, Li Shengyu, Lei Jiaqiang, et al. Sand flux estimation of 


different underlying surface, a case of Cele oasis-desert ecotone[J]. 
Arid Land Geography, 2017, 40(3): 533-540. ] 

李 汉 林 , 何 清 , 金 莉莉 . 塔克拉玛干 沙漠 腹地 和 北 缘 典 型 天 气 近 
地 层 风 速 廊 线 特征 转 . 干旱 气象 , 2020, 38(6): 965-978. [Li 


Hanlin, He Qing, Jin Lili. Wind profile characteristic of near-sur- 


face layer under typical weather conditions in the hinterland and 


202210.00036v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


northern margin of Taklimakan Desert[J]. Journal of Arid Meteorol- Technology, 2021, 23(4): 137-144. ] 
ogy, 2020, 38(6): 965-978. | [34] 何 清 , 胡 文 峰 , 杨 兴 华 , 等 . 巴 丹 吉 林 沙 漠 拐 子 湖 地 区 贴 地 层 风 
[32] MHAE, 张 克 存 , 安 志 山 , 等 . 敦 得 沙漠 绿洲 和 戈壁 地 表 风 动力 环 速 廊 线 和 风沙 流 结构 特征 团 . 干旱 区 研究 , 2012, 29(3): 517-523. 
境 特征 同步 对 比 加 . 中 国 沙漠 , 2020, 40(4): 113-119. [Hu Fei, [He Qing, Hu Wenfeng, Yang Xinghua, et al. Research on wind pro- 
Zhang Kecun, An Zhishan, et al. Comparison of wind dynamic en- file and sand drift structure in Guaizi Lake Region in the Badain 
vironment among desert, oasis and Gobi[J]. Journal of Desert Re- Jaran Desert[J]. Arid Zone Research, 2012, 29(3): 517-523. ] 
search, 2020, 40(4): 113-119. ] [35] 张 正 假 , 董 治 宝 , 赵 爱 国 , 等 . 沙漠 地 区 风沙 活动 特征 一 一 以 中 
[33] 贾 瑞 庭 , RINA, 党 晓 宏 , 等 . 库 布 齐 沙 漠 200 MWp 光伏 阵列 的 国 科 学 院 风沙 观测 场 为 例 困 .干旱 区 研究 , 2007, 24(4): 550- 
截流 阻 沙 效应 研究 趾 . 中 国 农业 科技 导报 , 2021, 23(4): 137- 555. [Zhang Zhengcai, Dong Zhibao, Zhao Aiguo, et al. Features 
144. [Jia Ruiting, Meng Zhongju, Dang Xiaohong, et al. Study on of sand drift movement in desert: A case study at sand drift obser- 
the intercepting and blocking sand effect of 200 MWp photovolta- vation station of Chinese Academy of Sciences[JJ. Arid Zone Re- 
ic array in Kubuqi Desert[J]. Journal of Agricultural Science and search, 2007, 24(4): 550-555. ] 


Aerodynamic parameters of typical underlying surfaces in an aeolian region in 


the middle reaches of the Yarlung Zangbo River 
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Abstract: Aerodynamic parameters are an important research focus in the study of aeolian areas. Some of the 
most severe aeolian damage occurs along the middle reaches of the Yarlung Zangbo River in Tibet. The 
floodplains and shifting sand dunes are the main surface types and source areas of wind-sand degradation in this 
region. The study of aerodynamic parameters is conducive to the assessment of aeolian damage and the 
forecasting of sandstorms. For this reason, we examined aeolian processes in three surface types: river floodplain 
vegetation area (vegetation coverage > 2096), disturbed river floodplain vegetation area (vegetation coverage < 
20%), and sand dune area (without any vegetation). The wind speed profile method was used to calculate the 
frictional wind speed (u:) and aerodynamic roughness (zo). The results showed the following: (1) The wind speed 
reduction amplitude near the ground layer was in the following order: vegetation area > disturbed vegetation 
area > sand dune area, and the wind speed reduction near the ground layer in the vegetation area was nearly twice 
that of the sand dune. (2) Although the near-surface wind speed profiles of different underlying surfaces can be 
expressed as log-linear functions (u; = bln(z) + a), the fitting parameters are different. (3) The averaged values of 
u» in the vegetation area, disturbed vegetation area, and sand dunes are 0.73, 0.59, and 0.36 m.s ', respectively, and 
the averaged values of zo are 7.1 x 107, 3.7 x 107, and 4.7 x 10* m, respectively. The u» and zo of the vegetation area 
are approximately 2 and 150 times those of the sand dunes, respectively. (4) The u» of the vegetation area, disturbed 
vegetation area, and sand dunes all increased with increasing averaged wind speed, and the increase in amplitude is 
in the following order: vegetation area > disturbed vegetation area > sand dunes. (5) The responses of zo to the 
average wind speed in the vegetation zone and the sand dunes are clearly different. The zo of the vegetation area 
increased slowly with the increase in wind speed. The value of z in the sand dunes area has a decreasing trend at low 
wind speeds but begins to increase when the wind speed increases to a threshold value of 6.5 m-s”. 


Keywords: Yarlung Zangbo River; wind speed profile; friction velocity; aerodynamics; roughness length 


